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The aim of this paper is to answer the following two questions: (1) Given cosmological observations 
of the expansion history and linear perturbations in a range of redshifts and scales as precise as is 
required, which of the properties of dark energy could actually be reconstructed without imposing 
any parameterization? (2) Are these observables sufficient to rule out not just a particular dark 
energy model, but the entire general class of viable models comprising a single scalar field? 

This paper bears both good and bad news. On one hand, we find that the goal of reconstructing 
dark energy models is fundamentally limited by the unobservability of the present values of the 
matter density flmo, the perturbation normalization as as well as the present matter power spectrum. 
On the other, we find that, under certain conditions, cosmological observations can nonetheless rule 
out the entire class of the most general single scalar-field models, i.e. those based on the Horndeski 
Lagrangian. 



I. INTRODUCTION 

Research in dark energy (DE) cosmology is gener- 
ally devoted to building viable models and to constrain- 
ing them from observations (see, for instance, the re- 
views [1-3]). The models are usually characterized by 
a small number of properties at background and linear- 
perturbation level (equations of state, speeds of sound, 
masses, coupling strengths, etc.) vsfhich then have some 
particular effect on the phenomenology (evolution of the 
scale factor and perturbations of matter and DE itself). 
In this paper, we employ a different approach, aiming to 
answer the following questions: 

1. Assuming the ideal case of cosmological observa- 
tions of the expansion history and linear perturba- 
tions in a range of redshifts and scales which are 
as precise as is required, which physical properties 
(such as e.g. the Hubble rate H{z), the perturba- 
tion normalization as, the perturbation growth rate 
/, etc.) could actually be reconstructed if we were 
to refrain from any parameterization of dark en- 
ergy? 

2. Can we use these observable quantities to rule out 
not just some particular cosmological model but 
the entire class of viable single scalar-field models? 

We are of course not the first to attempt to study dark en- 
ergy cosmologies in a model-independent way. In Ref. [4] 
the author argues for a reconstruction method which does 
not rely even on Einstein equations. This method can in 
principle directly measure the space-time curvature but 
cannot test a modification of gravity. On background 
level, the power of observations to constrain the expan- 
sion history of the universe by assuming A-cold dark mat- 
ter (ACDM) as a null test was studied in Refs [5, 6], while 
a principle-component analysis of the equation of state 



constraints was developed in Ref. [7]. Recently, even the 
validity of the assumption of the Copernican principle 
was tested [8]. On the level of linear perturbations, pa- 
rameterizations are usually used to limit the freedom in 
the model-independent description of growth of structure 
[9-13], although principle components analysis has also 
been employed [14]. Another way to limit the freedom 
in the fully general description is to exploit the struc- 
ture which any general-rclativity-like theory of dark en- 
ergy must obey: such approaches were discussed in [IS- 
IS]. All of the above approaches, when contrasted with 
data, require parameterizations in order to break degen- 
eracies, but simultaneously introduce parameterization- 
dependent biases. 

The difference of our approach is that, given the min- 
imum of assumptions, we first elucidate the observables 
that measurements can in principle provide without the 
assumption of any dark energy model in particular. It is 
only then that we use these model-independent observ- 
ables to construct tests which might eliminate or confirm 
particular models. Our approach is closest in spirit to 
Ref. [19], the results of which we extend. 

In this paper, we completely ignore the practical prob- 
lems and limitations of the observations and assume that 
good-enough statistics with sufficiently small systematic 
errors can be achieved in the range of redshifts and scales 
discussed here. By exploring this idealized case we try to 
discover the fundamental limits to which observations in 
a dark energy cosmology are subject. 

We adopt metric signature ( — h ++), while a comma 
denotes a partial derivative. We interchangeably use co- 
ordinate time t, scale factor a, e-foldings iV = Ina and 
the redshift z as time variables. Overdots denote deriva- 
tives with respect to t, primes with respect to N. The 
subscript denotes the present time. We also make use 
of the notation of [20, 21], to which we refer for a thor- 
ough study of linear perturbations in the context of the 
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Horndcski theories. 



II. ASSUMPTIONS 

In keeping with the spirit of gencrahty, we first wish 
to make the minimum of assumptions on the geometry 
and matter content of our Universe that wiU aUow us to 
interpret observations at aU. In the following we assume 
only that: 

(a) The geometry of the Universe is well described by 
small linear perturbations living in an FLRW met- 
ric with scale factor a{t). We will not consider pos- 
sible observations of rotational perturbation modes 
nor of gravitational waves, as these are irrelevant 
for structure formation in late-time cosmology. 

(6) The matter content (i.e. dark matter and baryonic 
matter) is pressureless or evolves in a known way. 

(c) The relation between the galaxy distribution and 
the matter distribution at linear scales can be mod- 
eled as (5gai = b{k, a)5ni, where 6(fc, a) is the poten- 
tially scale- and time-dependent linear bias, while 
at the same time there is no bias between the ve- 
locities of galaxies and matter. This implies that 
both the baryonic and dark matter respond in the 
same way to the gravitational potentials and that 
the statistical velocity bias due to galaxies sampling 
preferentially over-dense regions [22] is negligible on 
the scales of interest. 

(d) The late-time universe is effectively described by 
the action 




(setting SttGn = 1) which includes the Einstein- 
Hilbert term for the metric g^i, and the La- 
grangian £,„ describing pressureless matter fluids, 
both baryons and dark matter, between which we 
will not differentiate here. Any other terms are as- 
cribed to the DE Lagrangian , which represents 
some consistent theory potentially depending on 
extra degrees of freedom or g^i, (i.e. modifications 
of gravity).^ We will neglect the radiation com- 
ponent because all the observations arc assumed 
to be performed well after decoupling. In non- 
minimally coupled models, the Lagrangian de- 
pends on a different metric, related to g^j/ through 
some transformation. Here we assume, however, 
that we have already reformulated the action so 
that matter moves on the geodesies of 5^,^. 



^ A consistent theory is understood here to be a theory free of ghost 
and other catastrophic instabilities that can in general occur in 
generalized gravity and dark energy models. 



We employ the above minimal framework to address 
question (1) by considering the background observables 
in section III and those arising from linear perturbations 
in section IV. To answer question (2) we need another 
crucial assumption, concerning the degrees of freedom in 
the dark energy Lagrangian: 

(e) The Lagrangian Cx, which describes dark energy, 
is any one of the Lagrangians describing a single 
scalar field governed by second-order equations of 
motion. We call this scalar field dark energy, but 
we do not necessarily require it to be driving the 
current acceleration. For example, it could be that 
the late time acceleration is effectively driven by a 
cosmological constant, but in the presence of this 
additional degree of freedom. The assumption of 
a scalar field ensures that there are no gross viola- 
tions of isotropy. The limitation to second order is 
a necessary condition to ensure that the model is 
not subject to instabilities (see e.g. [23]).^ We will 
therefore assume that the dark energy is governed 
by the most general Lagrangian which fulfills these 
requirements: + Cx will form the Horndeski La- 
grangian (HL [25, 26]). We dedicate section V to 
this system. 

III. BACKGROUND OBSERVABLES 

From assumptions (a)-(c), by varying the action eq. (1) 
with respect to the metric, we obtain a Fricdniann equa- 
tion that can be written as 

where iJg is the present value of the Hubble parame- 
ter, f2fco the present curvature density parameter and pm 
is the matter energy density. From assumption (6), 
evolves as , and px is the energy density of the terms 
coming from Cx- 

Observations of the cosmic expansion are essentially 
estimations of distances D{z) (i.e. luminosity or angular- 
diameter distances) or directly H{z) (e.g. using mea- 
surements of longitudinal baryon acoustic oscillations, or 
real-time redshift-drift observations [27]) based on the 
existence of standard candles, rods or clocks. More ex- 
actly, standard candles or rods measure H{z) up to a 
multiplicative constant, related to the unknown absolute 
measure of the source luminosity or proper length. For 
instance, the flux of supernovae la (SNIa) with absolute 
luminosity L are known only up to the constant LH^; 
only ratios of fluxes at different redshifts are independent 



^ This class of Lagrangians includes such theories as f{R) grav- 
ity, despite their naively fourth-order equations of motion. This 
is because we can always introduce the a priori hidden scalar 
explicitly through a Legendre transformation [24]. 
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of the absolute normalization. The same is true of baryon 
acoustic oscillations: they measure only the ratio of the 
sound horizon at last scattering and the Hubble radius 
H~^{z) . We can therefore say that, without additional 
assumptions, background cosmological observations esti- 
mate D{z) up to an overall constant as well as the di- 
mensionless Hubble function E{z) = H{z)/Ho. Notice, 
however, that real-time redshift-drift observations can es- 
timate the absolute value of H{z), while local measure- 
ments of the expansion measure Hq. 

Combining D(z) with H(z), we can also estimate the 
present curvature parameter ^ko- We can therefore de- 
termine the evolution of the combined matter and dark 
energy content, 1 — 51^, at all times. If we assume that 
there are only two components of the cosmic fluid then 
we have only one free parameter, flmo- In fact, we can 
write 



(3) 



Therefore, we conclude that from background observ- 
ables we can reconstruct both and but only up 
to r^mo [28] , since one can compensate for any change of 
rimo with a modification of the DE model. Of course, if 
we parameterize the evolution of fl^ with a simple equa- 
tion of state, we can break the degeneracy with finiOi is 
usually done in analyses of SNIa data, but that is exactly 
what we are trying to avoid in this work. 

The same result is valid if instead of pure pressure- 
less matter one includes further components (e.g. mas- 
sive neutrinos) that evolve with an effective equation of 
state Winl-z), provided Wmiz) can be inferred from other 
observations (e.g. knowledge of the neutrino masses). 



IV. LINEAR PERTURBATION OBSERVABLES 

The linear perturbation observables are the correla- 
tions of positions, velocities and shapes (ellipticitics) 
of sources (i.e. galaxies, Lyman-a lines, clusters, back- 
ground radiation) in angular separation and redshift. 
Given knowledge of D{z), these can be converted to the 
more usual dependence on wavenumber k and redshift. 

Let us first discuss the clustering of matter. We denote 
the root mean square of the correlation of galaxy num- 

1 /2 

ber counts in Fourier space as (Sgai (i.e. Sg^i = P^^i (fc, z) 
where Pgai is the galaxy power spectrum). We de- 
fine from now on the wavenumber k to be the physi- 
cal wavenumber expressed in the units of the cosmolog- 
ical horizon, i.e. h — ^phys 

/ aH (fc is independent of Hq 
if fcphys is measured in h/Mpc). This means that k is 
time-dependent. We observe galaxies, not matter per- 
turbations, so as anticipated we need to introduce a bias 
function b{k,z) such that (5gai = bSm- Without the as- 
sumption of a particular model, DE perturbations are 
unknown. In many models they are not at all small com- 
pared to matter perturbations. For this reason, we will 



define the total density perturbation St = ^^m'^m + ^x^x 
and introduce the bias B of galaxies with respect to it. 



(5gal = B5t = BZVtraK 



(4) 



where Z{k,a) = 1 + n^Sx/ {nraSra) is a function of 
space and time that depends on the clustering of the 
x-component. Then we have that b = BZfl^. 

Let us denote the initial total density perturbation 
spectrum at decoupling as 5t^j^{k) and as Gt(k,z) the 
scale-dependent growth function of the linear total den- 
sity perturbations, normalized to unity at present. If 
the galaxies move with the same velocity field as mat- 
ter, the galaxy velocity divergence 6'gai in the sub-Hubble 
regime is related to the matter density perturbation as 
6'gai = (^m = = ^f^ra, by the Continuity equation 

for matter. We then obtain 0gai = —{f/b)5gah where 
/ = G' /G is the linear matter growth rate and G(fc, z) 
is the growth function for matter perturbations, both of 
which are scale-dependent for general DE models. This 
velocity field generates redshift distortions as a function 
of the direction cosine fi ~ {k ■ £)/k where £ is the unit 
line-of-sight vector. The observable (5gai can therefore be 
expressed as [29] 



(5gai(fc,z,^) = GtBcrs,t 



1 + p 



StAk), (5) 



where crg t is the present normalization of the to- 
tal density spectrum. Now we can write GtB = 
Zilin(5nii?/(Zr2in(5m)o = Gfe/(Znni)o, SO wc havc the 
almost-standard expression 



(5gai(A:,z,^) = Gbas 1 



/ 2 
-b^ 



'5t,o(fc) 



(6) 



In this expression we set, using the 8 Mpc/h spherical 
Fourier space window function W{k), 



CT8 



{zn 



Zr. 



m)0 



1/2 

Z'^P,^oW'{k,Rs)dk] ,(7) 



which is equivalent to the usual normalization ag, if Z 
depends weakly on k. 

To discuss weak lensing, we introduce the standard 
perturbed metric in longitudinal gauge ds^ = —(1-1- 
2^')dt^ -I- a^(l + 2$)(5yda;Mx-''. Since later we wiU con- 
fine ourselves to a scalar dark energy, as part of assump- 
tion (a), we assume that the DE does not excite vec- 
tor and tensor modes so that only scalar modes need to 
be included. It is helpful to introduce the function Y, 
the effective gravitational constant for matter, and the 
anisotropic stress r] defined as (see e.g. [12, 21]) 



Y{k,a) 



2fc2* 



ri{k,a) 



(8) 



Both rj and Y are unity on sub-Hubble scales if Cx is a 
constant (i.e. the DE is a cosmological constant) and the 
matter is a perfect fluid. In the linear regime, the lensing 
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effect is proportional to the Icnsing potential, which itself 
is driven by the density perturbations (see e.g. [30]). In 
general, this relation can be written as 

fc2$i,„3 = fc2(^ _ $) ^ _'^Y{1 + 7j)n,^S,-^ (9) 

where we have defined the "modified lensing" function 
z) = Y{1 + T]). The ellipticity correlation is an in- 
tegral function of within a window function that 
depends on the survey geometry (see e.g. [19]). Assum- 
ing good-enough knowledge of the galaxy distribution one 
can differentiate the correlation integral and obtain the 
quantity 

a{k, z) = ^(fc4p<j,^_)i/2 ^ 1 j^^^pSG^gj . (10) 

Then from Sgi^^k, z, fj.) (with e.g. /i = 0, 1) and (j{k,z) 
one can measure the three quantities R, L defined as 

A = GbasStfi , R = Gfas5t,o , (11) 
L = n^o^GasSt.o ■ (12) 

The quantities that connect the observations to theory 
(i.e. to the Lagrangian Cx) are Q^i /, so it would 
be optimal to estimate them directly from observations. 
Now, the cosmic microwave background anisotropy al- 
lows one to measure, at least in principle, the initial po- 
tential ^Pin through the Sachs- Wolfe effect. It is, however, 
impossible to derive from this information the present 
power spectrum 6tfi since it also depends on a scale- and 
time-dependent transfer function. Absent a model for 
DE, this transfer function is unknown and since it acts 
to process the total perturbation spectrum, changing its 
fc-dependence, it also makes St^ an unknown without fur- 
ther assumptions. This argument shows that the only 
(5t,o-independent quantities directly measurable from lin- 
ear cosmological observations are ratios of A, R, L, and 
their A^-derivatives, i.e. 

Pi = R/A = (13) 
P2 = L/R = a„oS//, (14) 
Pa = R'/R = / + /'//. (15) 

All other possible (5t,o-independent ratios, such as A' /A, 
L' /L or R' / L or higher-order A^-derivatives, can be ob- 
tained as combinations of P1-3 and their derivatives, for 
instance L' /L = P2/P2 + Pa- Other linear-perturbation 
probes, such as integrated Sachs- Wolfe, cross-correlations 
or 21-cm flux measurements, add statistics and might ex- 
tend the observational range but do not break the fun- 
damental degeneracy. 

The quantity Pi is a well-known observable quantity, 
often denoted /3 [30]; since it involves the bias function 
&, related in an unknown way to the model of dark en- 
ergy, we will not consider it any longer in this paper. 
The quantity P2 has already been introduced in [19] as 



Eq as a test of modified gravity, but the fact that flaio 
is not an observable was not discussed there. The quan- 
tity R contains the term Gfas, also denoted as f<Js{z) 
in the literature [31]. This term is often considered to be 
a directly observable quantity, but, as we have argued, 
this is only true if one assumes a model for DE, or at 
least a parametrized form of St.o] otherwise, the model- 
independent observable combination is P3 — R'/R. It is 
important to realize that even a perfect knowledge of 
P3 does not imply knowledge of / since the equation 
/'// + / = Paif^jz) cannot be solved without the un- 
known /c-dependent initial condition for /. Finally, no- 
tice that we did not need to assume Gaussian fluctuations 
nor isotropy of the power spectrum. 

Measurements of galaxy peculiar velocities and their 
time derivative directly estimate through the Eu- 
Icr equation, which would give the quantity V — 
riijioYGasStfi ■ Then one can form the observable, 
L/V = 1 + 1], which measures the anisotropic stress. This 
new observable is not independent since it can be written 
in terms of E,P2,P3 (see Eq. (21)). Moreover, the esti- 
mation of V requires a delicate subtraction of the peculiar 
redshift from the cosmological redshift by using distance 
indicators such as Cepheids and therefore a number of ad- 
ditional assumptions on the source physics. No current or 
foreseeable method to estimate the peculiar velocity field 
(let alone its derivative) has been shown to be reliable 
beyond a few hundred megaparsecs (see e.g. Ref. [32]), 
so we will not pursue this possibility any further in this 
paper. 

Our first result is that linear cosmological observations 
can at best determine only E = H/Hq as function of time 
(as well as Vlko but not Qrno) a-nd the observable combi- 
nations Pi_3 as functions of time and space, within the 
range of the observations themselves. To achieve this, we 
need to combine galaxy clustering and weak lensing in the 
same redshift range. It is possible that one can determine 
other combinations from non-linear effects, but this will 
certainly introduce new uncertainties (e.g. a non-linear 
bias). 



V. THE HORNDESKI LAGRANGIAN 

It is now time to use also assumption (e) regarding 
the explicit form of the dark energy Lagrangian. In the 
choice of the action (1), we have assumed that all mat- 
ter components (i.e. dark matter and baryons) feel the 
same gravitational force and propagate on geodesies of 
the metric g^i,. In addition, we now explicitly assume 
that the DE is modeled by a single scalar field cj) described 
by the Horndeski Lagrangian (HL). The HL is defined as 
the sum of four terms £2 to £5 that are fully specified 
by a non-canonical kinetic term K{(j},X) and three in 
principle arbitrary coupling functions 03^4,5(0, A), where 
X = — (7^iy(/)'^(/)''^/2 is the canonical kinetic term. 
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£2 =K{<P,X), 



£5 =G5(<^,X)G^,V^V''0- 



G 



(16) 



5,X 



6 

2(V^V^v 



(□0) - 
3" 



The Horndeski Lagrangian is the most general La- 
grangian for a single scalar which gives second-order 
equations of motion for both the scalar and the met- 
ric on an arbitrary background. This is a necessary, but 
not sufhcient, condition for the absence of ghosts. In 
general, the equation of motion for the scalar will couple 
it to the matter energy density. The metric potentials 
$ and 4" are as usual determined by the Poisson and 
anisotropy equations, which are constraints, and there- 
fore do not have independent dynamics. We note that 
the generalization to the case of multiple scalar fields has 
been discussed in Ref. [33]. 

In what follows, we will assume that the so-called 
quasi-static limit is valid for the evolution of pertur- 
bations. This implies that we are observing scales 
significantly inside the cosmological horizon, k = 
kphys/ (aH) 3> 1, and inside the Jeans length of the 
scalar, Cgk ^ 1, such that the terms containing k dom- 
inate over the time-derivative terms. The sound speed 
Cs is a particular function of the HL functions K, G3_5 
evaluated at the background level and we have presented 
it in Eq. (A4) [20]. In this quasi-static limit for a model 
belonging to the HL, one obtains [21] 



1] = /l2 



1 + k^h. 



(17) 



for suitably defined functions /11-5. In this limit, one also 
has Z = Yij. 

The functions /ii_5 express the modification of gravity 
induced by the HL. In real space they induce a time- 
dependent Yukawa correction to the Newtonian poten- 
tial. They are all combinations of HL functions K, G3,4_,5 
and their derivatives with respect to (j) and X, all evalu- 
ated on the background and arc therefore time- but not 
/c-dependent, 

hi EE hi{z) EE h,{4>,X). 

The explicit expressions for the functions hi are very 
complicated and not particularly illuminating; we have 
nonetheless presented them in appendix A. For ACDM 
one has simply hi_2 = 1 and /i3,4.5 = 0. If the two 



^ The constraints on the HL arising from stability considerations 
are derived in Ref. [20] . We have also presented them in appendix 



gravity-coupling functions in the HL, 0^,0^, are con- 
stant, i.e. the effective Planck mass is constant, then 
r] — 1; if moreover G3 depends only on (i.e. k-essence), 
then also Y = 1 and there are no modified-gravity effects 
at all in this quasi-static limit. 

It is worth noting that one could have arrived at the 
form of Eqs (17) given just our assumptions of second- 
order equations of motion, the symmetries of the FRW 
background and quasi-staticity. 

In the same quasi-static limit, from the matter conser- 
vation equation, we obtain 



S'i + 2 + 



hi 



1 



(18) 



On the other hand, we can write for the weak-lensing 
function S, 



s = y(i + 77) = 



1 + k'^hj 



(20) 



where we have introduced two auxiliary functions, /ig = 
hi[\ + /12) and hj = {h^ + /i4ft.2)/(l + ^2)- As an aside, 
one can show that if G3 , G4 depend only on (p and G5 = 
const, i.e. we are dealing with a k- essence theory non- 
minimally coupled to gravity, then /12 = 1 and /17 = ^13 
so that S becomes independent of fc. In this limit, the 
gravitational potential felt by photons is not distorted 
(we have discussed this model in detail in Ref. [34]). 

From the observables Pi, E we can construct a model- 
independent relation measuring the anisotropic stress rj 
as follows. From ^2,^3 we can obtain / ~ ri,„oS/-P2 
and /' = P3f}moi;/P2 - (l^moS/P2)^- Inserting this in 
Eq. (19) and employing Eq. (20), after a little algebra we 
obtain a simple relation. 



3P2(l + 2) = 



2E^ R 



ML' 

E . 



- 1 = ^7 = ^2 



It is important to stress that the l.h.s. of Eq. (21) is a 
function of model-independent observables, and thus a 
model-independent measurement of ry, valid under our 
assumptions (a)-(c?), but not requiring (e). The form of 
the last term is determined by the QS limit of the Horn- 
deski Lagrangian. Given the above, we can exclude all 
dark energy models described by a single scalar field in 
the QS limit by showing that the anisotropic stress mea- 
sured from the observation data does not follow the par- 
ticular scale dependence mandated by Eq. (21). Equa- 
tion (21) must be valid in fact at all times and scales 
where the quasi-static limit is valid. At any given epoch 
z* , this equation involves the three unknowns /i2,^4,^5 
all evaluated at 2*. If at this epoch we observe E and 
P2, P3 at more than three fc-modes, we can form an over- 
constrained system. If for any z* this system has no 
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solution then the observations are inconsistent with the 
quasi-static Hmit of HL. Equivalently, from Eq. (21) one 
can obtain a consistency relation that depends only on 

observable quantities. Defining g{z,k) = ^"^^°2^ ; one 
has in fact 

2.9^'\g('^-3(.g(2))2 = o, (22) 

where the n-th derivative of g with respect to k . 

If this condition fails at any one redshift, the DE is not 
described by the HL in the linear quasi-static limit. This 
is the second main result of this paper. Needless to say, a 
cosmological constant satisfies this consistency relation. 

On the other hand, if there are consistent solutions 
then we obtain an indication in favor of the HL and 
also direct constraints on it. For instance, if observa- 
tionally we find that P2,^3 do not depend on k, then 
from Eq. (21) we see that the condition ^,4 ~ /15 must be 
satisfied. 

If the consistency relation is not satisfied, the only pos- 
sible way out of our conclusion is that the conditions 
for the linear quasi-static limit that we employed to de- 
rive Eq. (17) are not satisfied. This can occur if the 
rate of change of the functions /ii-5 is very large, e.g. if 
h'j/hj w Cg/c^ ^ 1 for some j. However if the field cj) 
drives the current accelerated expansion (this is indeed 
an additional assumption) we expect it to be slow rolling 
on time scales of the order of hj/hj ^ H, i.e. h'j/hj w 1. 
If the sound speed squared is of order unity, then 
on typical astrophysical scales of 100 Mpc/ft. one has 
k^ « 10"^, so the quasi-static limit should be very well 
satisfied. However if Cg is less than say 10~^ then the 
simple form of Eq. (17) is no longer valid. One might 
then expect oscillating terms in the Y, rj equations; it is 
possible that this behavior could be probed, and possi- 
bly rejected, by a similar method we are discussing here 
but a full analysis of this "cold dark energy" scenario 
would be required. Another potential difficulty is the 
fact that these more general scalar field theories contain 
non-linearities in principle independent of those in the 
matter perturbations. It may prove difficult to deter- 
mine on which scales the linear approximation for the 
dark energy is valid, if at all. 

From P2 and P3 one can build other consistency 
equations, e.g. by differentiating P2 or the combination 
P2/P2 + -P3 = S'/S -I- / with respect to N and again em- 
ploying Eqs (19), (20). These relations however require 
derivatives of the observables P2 , P3 and will introduce 
derivatives of the /11-5 functions, so appear to be less 



useful than Eq. (21). 

We observe also that the propagation speed ct of 
gravitational waves is a function of the HL coupling 
fimctions, see Eq. (A3) [35]. A detection of a source both 
in gravitational and electromagnetic waves could allow 
for a measurement of Ct and therefore new independent 
constraints on the HL [36] . 

VI. CONCLUSIONS 

We have shown that cosmological linear observations 
can measure only flko, E = H/Hq and the combinations 
Pi = f/b, P2 = f^moS// and P3 = / + /'//■ Parameters 
such as riniOj or functions such as Gfas are not di- 
rectly model-independent measurable quantities via lin- 
ear cosmological observations alone. This limits in a fun- 
damental way the knowledge of, among others, the evo- 
lution of the DE density parameter Q^, its equation of 
state, or the matter growth rate /. From _B, P2.3 one can 
form consistency relations in terms of the HL functions. 
The simplest one is Eq. (21) or (22), expressed purely in 
terms of model- independent observables. If observations 
indicate a violation of a consistency relation, then the 
DE is not described by the HL in the quasi-static limit. 
Conversely, finding the predicted fc-behavior would be a 
major confirmation of the scalar field picture of dark en- 
ergy 

Non-linear effects will bring both more information 
and more unknowns into the picture so it is not clear 
how much, if at all, they would improve the task of 
reconstructing or rejecting the HL. The limitations of 
real-world observations, completely neglected here, are 
of course in practice the major hurdle on the path to this 
goal. 
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Appendix A: Details of Horndeski Properties 



This appendix concerns the properties of the scalar-field theories described by the Horndeski Lagrangian given by 
the combination of terms presented in Eqs (16). The subscripts , (/> and ,X denote derivation with respect to that 
variable. On a flat FRW background, the energy density and pressure are given by 
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p., =3H^{1 ~wi) + 2XK,x -K~ 2XG3,^+ (Al) 

+ 6^H {XG3,x - G4,0 - 2XGi^^x) + 

+ 12H {X (6*4, A' + '^XGi^xx) — G5,0 — XG^^c^x) + 

+ A^XH^ {G5,x + XGz ,xx) , 
P^=- (3ff2 + 2H){l-wi)+K- 2XG3,0 + 4XG4,00+ 

+ 2<pHwi^^ - iX^H^G^^^x + 2<pXH^G5,x + t (^2 - 2Hwi) , 

where, given a slight rearrangement of the resuhs in Refs [20, 21], we define four functions Wi as 

w;i =1 + 2 (G4 - 2XG^,x + XGs,^ - ^XHG^.x) , (A2) 
W2=-2^ {XG^,x - G4,0 - 2XG4,0x) + 

+ 2H {wi - AX (G4,x + 2XGi,xx - Gg^^ - XG^,^x)) - 

— 2^XH^ {3G5.X + 2XG5 ,xx) , 

^3 =3X {K^x + 2XK ,xx — 2G3_0 — 2XG3^ct>x ) + 18<j>XH {2G3,x + XG3 ,xx) — 

— 18(j>H (G4,0 + 5XG4^cj)X + 2X G4^^xx) — 

— 18-ff" (1 + G4 — 7^G4.x — 16^'^G4,xx — 4Jf '^G4^xxx) — 

— 18Xi?^^6G5_0 + 9XG5,0x + 2Ar^G5_0xx^ + 

+ 60Xi/3(^15G5,x + 13XG5,xx + 2X^G5,xxx) , 
W4 =1 + 2 (G4 ' XG5,0 - XG5^x4>) ■ 

All of the dynamics of hnear perturbations are fuUy determined by the above four functions. In particular, the speed 
of propagation of gravitational waves, ct, and the normalization of the kinetic term of these tensor perturbations, 
Qt, is given by 

4 = !^^>0, Qt=^>0, (A3) 
wi 4 

with positivity required by stability. From the above, is can be seen that wi has the meaning of the normalization 
of the tensor perturbations, i.e. it is the effective Planck mass squared. The corresponding quantities for the scalar 
degree of freedom, the sound speed of dark energy, Cg, and the normalization of the kinetic energy for perturbations, 
Qs, in the presence of dust with energy density pm, are 

2 _ 3 {2wlw2H - W|W4 + AW1W2W1 - 2wl{w2 + Pm)) , . , 



wi (AwiWs + 9m 



With above definitions in hand, we can define the five scale-independent functions /ii_5 which appeared in the result 
(17). All the observables for scalar perturbations in the quasi-static regime are determined by these five functions, 

(A5) 



hi = 


W4 


_Wl _ 








WI Wi 


^ W4 








2wlw2H 


— W2W4 + 


AW1W2WI 


- 2wl{w2 + Pm) 


h3 = 


2XAP 




2wl 




/14 = 


2wfH^ ~ 


W2W4H 4 


- 2wiWiH 


+ W2W1 - Wi(w2 + Pm) 




2XM^ 




Wl 






2wlH^ - 


■W2W4H 4 


- AwiWiH 


+ 2wi'^ - W4,{W2 + Pm) 


/15 = 


2X1VP 




W4 
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and where the effective mass squared, M^, can be expressed in terms of derivatives of the total pressure and total 
energy with respect to the scalar as 

One may wonder whether it is possible to invert Eqs (A5) in order to obtain pm as a function of the functions /ii_5. 
If this were possible, and if all the /ii_5 were observable, then one could measure fimo, contrary to our claim in the 
text. In a future paper we will discuss in detail the observability of the functions /ii_5 and we will show that in fact 
it is not possible to obtain ri,„o this way. Here it will suffice to notice that hi is completely degenerate with i7,„o (see 
Eq. (19)); it turns out that p,,! is proportional to 1/hi and therefore the fact that hi can only be measured up to fimo 
implies the same degeneracy in pm- 
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